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INTRODUCTION 

The  function  of  proteins  varies  widely  in  nature  as  a  result  of  evolution  and 
includes  proteins  which  are  toxic  to  cells.  The  protein  ricin  is  among  the  most  cytotoxic 
substances  known  [1,2].  Ricin  functions  by  inhibiting  the  synthesis  of  protein  in 
eukaryotic  cells  through  the  hydrolysis  of  an  N-glycosidic  bond  of  the  adenosine 
residue  in  28S  ribosomal  RNA  (rRNA)  [3.4].  Ricin,  containing  the  A-  and  B-chains,  is 
isolated  from  the  seeds  of  the  castor  plant,  Ricinus  communis.  The  A-chain  (RTA)  is 
composed  of  267  amino  acid  residues  and  contains  the  enzymatic  activity;  the  B-cham 
(RTB)  is  composed  of  262  residues  and  binds  to  galactose  residues  present  on 
various  cell  surface  glycoproteins  and  glycolipids.  Depending  on  the  source  of 
ribosomes  and  assay  conditions,  RTA  has  a  and  Km  in  the  range  of  300*1500  per 

min  and  0.1-1 .3  pM.  respectively  [5.6].  It  has  been  determined  recently  that  RTA  will 

catalyze  the  deadenylation  of  a  synthetic  GAGA  tetraloop  oligoribonucleotide  with  a 
helical  stem  of  three  base-pairs  [7].  The  nature  of  the  base-pairs  does  not  influence 
recognition  or  catalysis  by  RTA,  whereas  the  position  of  the  tetranucleotide  in  the  loop 
sequence  is  critical  [7], 

In  an  attempt  to  elucidate  the  structural  interactions  governing  the  binding  of 
ricin  with  ribosomes,  Monzingo  and  Robertus  [8]  have  resolved  the  three-dimensional 
crystal  structures  of  formycin  S'-monophosphate  (FMP)  (structure  I)  and  adenyl-3’,5'- 

guanosme  (ApG)  (structure  II)  bound  in  the  active  site  of  RTA.  Structures  were  refintd 
at  2.8  A  and  3.0  A  resolution,  respectively.  These  two  nucleotide  ligands  are  thought 

m ,."\c  cer!a  ~  e  ements  ot  the  rRNA  substrate,  m  particuiar  the  bma.ng  of  tne 
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adenine  ring.  The  bound  ligands  lie  perpendicular  to  a  long  connecting  a-helix  near 
Glu  177  and  Arg  180  with  the  adenine  and  pyrozolopyrimidine  rings  stacked  between 
a  pair  of  tyrosines,  one  from  the  N-terminus  of  an  adjacent  p-strand  and  the  second 

from  an  a-helix.  The  A-chain  makes  additional  contacts  through  highly  polar  and 
charged  residues  located  within  joining  a-helices  linking  domains  two  and  throe. 
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To  date,  no  small  molecule  inhibitors  are  known  for  ricm  or  any  ot  the 
homologous  plant  or  bacterial  toxins.  Schramm  and  coworkers  [9.10]  have  shown 
FMP  to  be  a  powerful  transition-state  inhibitor  of  AMP  nucleosidase  from  Azotobacter 
vmelandii.  an  enzyme  that  hydrolyses  the  adenine  base  from  AMP.  The  tight  binding 
and  inhibitory  activity  of  FMP  has  been  attributed  to  a  binding  mode  of  locking  the 
pyrazolopyrimidine  ring  in  a  syn  configuration  about  its  C-glycosyl  bond  [11].  Ricm  A- 
chain  recognizes  and  binds  FMP  in  a  syn  configuration,  but  is  not  enzymatically 
inhibited,  and  in  fact  forms  a  weak  binding  complex  [8]. 

Analysis  by  X-ray  crystallography  suggests  that  ApG  might  bind  less  than  FMP 
[8].  Thus,  a  knowledge  of  the  structural  and  energetic  determinants  of  binding  ApG 
versus  FMP  to  the  A-chain  is  clearly  important  if  small  nucleotides  are  to  provide  a 
template  for  successful  inhibitor  design.  Moreover,  a  detailed  understanding  of  the 
conformational  requisites  that  binds  a  GAGA  tetraloop  to  the  active  site  of  ricin  A-chain 
is  likely  to  be  of  great  assistance  in  the  development  of  selective  inhibitors.  Towards 
this  goal,  molecular-dynamics  (MO)  simulation  studies  have  been  performed 
investigating  the  binding  of  FMP,  ApG.  and  the  hexanucleotide  CGAGAG  (data 
analysis  currently  in  progress)  with  ricin  A-chain.  In  addition,  complementary 
calculations  have  been  carried  out  aimed  at  exploring  various  structural  motifs  of  FMP 
which  might  increase  the  binding  affinity  for  the  active  site  as  well  as  the  da  novo 
design  of  small  ricin  inhibitors  employing  FMP  as  a  template. 

METHODS.  RESULTS  AND  DISCUSSION 

In  this  section  we  present  the  computational  methods  employed  and  the  results 
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obtained  for  the  molecular  modeling  studies  outlined  above.  Discussion  of  these 
relative  to  the  goals  of  the  research  is  presented. 

Molecular  Modeling  Methodology 

Molecular-dynamics  simulations  were  carried  out  using  the  program  DISCOVER  (12] 
with  force  fields  modeled  using  either  AMBER  (13)  or  CVFF  [14].  Atomic  charges  for 
FMP,  ApG.  and  CGAGAG  substrates  were  determined  by  either  electrostatic  potential 
fitting  to  a  single-point  wave  function  with  a  6-31 G*  basis  set  using  Gaussian  92  (15]  or 
from  the  AMBER  (CVFF)  nucleic  acid  library.  Parameters  for  the  protein  were  taken 
from  the  AMBER  or  CVFF  amino  acid  residue  library. 

The  substrate-bound  protein  structures  were  simulated  by  constructing  an 
active-site  region  composed  of  part  of  the  protein  and  substrate  inside  a  12-A 
spherical  boundary  as  derived  from  the  various  crystal  structures  (FMP  and  ApG)  and 
from  a  model  docking  CGAGAG  in  the  active-site  cavity.  Amino  acid  residues  lying  in 
the  outer  shed  were  fixed  at  their  initial  positions  and  should  not  introduce  significant 
phase-space  sampling  errors.  Active- site  regions  for  the  complexed  structures  were 
immersed  in  a  8-A  deformable  layer  of  water.  Water  was  represented  using  a  flexible 
SPC  potential  [16].  Simulations  of  the  various  RTA-substrate  active  regions  consisted 
of:  (1)  FMP  (or  analog)  plus  63  RTA  residues  (1026  protein  atoms),  two  sodium 
counterions,  and  362  water  molecules;  (2)  ApG  plus  83  residues  (1396  atoms),  one 
counterion,  and  544  water  molecules:  and  (3)  CGAGAG  plus  144  residues  (2350 
atoms),  six  counterions,  and  1231  solvent  molecules. 

Simulations  were  initiated  with  2000  cycles  of  minimization  using  a  steepest 
descent  algorithm  followed  by  30-ps  equilibration.  The  initial  atomic  velocites  were 
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assiqnea  from  a  Gaussian  distribution  corresponding  to  a  temperature  of  3C0  K. 
Nonbondei  mteractions  were  smoothed  to  zero  beyond  9.0  A  and  a  constant 
dielectric  (c  =  1 )  was  used.  All  hydrogens  were  treated  explicitly  employing  a  1.0  fs 

timestep  for  integrating  the  equations  of  motion.  Ensemble  averages  were  determined 
from  a  production  runs  of  100-300  ps  with  coordinates,  velocities,  and  energies  saved 
every  SO  timesteps  (0.05  ps)  for  further  analysis.  All  simulations  were  carried  out  on  a 
CRAY  Y-MP  located  at  the  National  Cancer  Institute.  Frederick  Cancer  Research  and 
Development  Center. 

Several  structural  motifs  constructed  from  the  original  FMP-bound  structure 
were  investigated  employing  free-energy  simulation  methods  [17].  For  estimating  the 
relative  binding  affinity  of  FMP  and  its  analogs  (denoted  as  FMP')  with  the  active  site  of 

RTA.  we  make  use  of  the  following  thermodynamic  cycle:  A(A A)  =  AA  intub’  -  A Ainh,b  = 
A At**  •  aAmw.  where  AA«**  and  A A*#*  are  the  Helmholtz  free  energy  of  binding  for 
substrates  FMP  and  FMP*.  respectively.  A  A**  represents  the  solvation  free-energy 
difference  between  FMP  and  FMP*.  and  A represents  the  binding  free-energy 
difference  between  FMP  and  FMP'  in  the  RTA  active  site.  Molecular-dynamics 
simulations  were  performed  for  a  series  of  small  discrete  transformations  of  FMP  into 
the  2*  and  2  -ammo  and  2-hydroxyl  analogs.  The  force  fields  were  modeled  using  the 

CVFF  calibrated  via  MNDO  [18]  calculations.  Simulations  of  the  substrate-bound 

e-ry—e  .  ac:  •.  a-'eg  s->  reth^dc’cgy  as  descr  ted  atc-.e  ^ 2 
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unbound  ligand  systems,  simulations  were  carried  out  by  placing  the  solute  from  the  X- 
ray  crystal  structure  in  a  cubic  box  of  length  27  A  containing  632  water  molecules  plus 
counterions.  Periodic  boundary  conditions  were  utilized  in  these  latter  set  of 
calculations  in  reducing  edge  effects. 

Free-energy  simulation  protocol  consisted  of  all  systems  simulated  at  ten  values 
of  the  advancement  coordinate  X  connecting  the  initial  (crystal  structure)  state  to  the 

mutated  state.  The  Xi  simulations  were  initiated  with  400  cycles  of  minimization 
using  a  steepest  descent  algorithm  plus  1000  cycles  of  minimization  via  a  conjugate 

gradient  algorithm  followed  by  10-ps  equilibration.  Simulations  for  X2,  X3 . X9  were 

carried  out  following  a  5-ps  equilibration  started  from  the  final  step  of  the  previous 
simulation.  The  Xio  simulations  were  initiated  with  a  10-ps  equilibration. 

The  de  novo  design  of  ricin  inhibitors  was  performed  using  a  computer  program 
called  LUDI  [19J.  As  described  in  the  literature,  LUDI  can  append  new  substituents 
onto  an  already  existing  ligand  [e.g.,  FMP)  in  such  a  way  that  hydrogen  bonds  are 
formed  with  the  protein  and  hydrophobic  pockets  are  filled  with  suitable  side  chains  of 
the  ligand.  Current  efforts  are  aimed  at  designing  new  substituents  for  the 
pyrazolopyrimidine  ring  given  the  fact  that  ricin  A-chain  recognizes  both  FMP  and 
ApG.  Energy  refinement  of  the  de  novo  designed  complexes  were  carried  out  via 
energy  minimization.  Poor  van  der  Waals  contacts  with  the  protein  and  poor  internal 
geometries  from  linking  multiple  fragments  are  removed.  Measures  of  the  fit  and 
binding  of  ligands  were  evaluated  through  the  total  nonbonded  interaction  erarg:°s 
with  the  receptor. 
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Results  and  discussion 


Summarized  in  Table  1  are  the  averaged  intramolecular  and  nonbonded  interaction 
energies  tor  the  RTA-FMP  complex  computed  over  200  ps  for  the  simulation  model 
(AMBER  force  fields)  and  the  corresponding  values  for  the  X-ray  structure.  The  energy 
contribution  that  shows  significant  deviation  is  the  nonbonded  electrostatic  interactions 
between  RTA  and  FMP.  Superposition  of  the  X-ray  structure  and  simulation  structure 
(shown  in  Fig.  1)  illustrates  key  structural  similarities  as  well  as  noticeable  differences. 
In  accord  with  the  crystal  structure,  the  simulation  structure  shows  the  formycin  ring 
conformationally  stacked  between  the  rings  of  tyrosines  80  and  123.  These  jo 
residues  make  strong  packing  interactions  with  the  substrate  stabilizing  the  binding 
complex.  The  hydroxyl  group  of  Tyr  80  further  stabilizes  the  binding  of  FMP  by  making 
several  distant  electrostatic  interactions  with  the  pyrazolopyrimidine  moiety  and  the 
ribose. 

Both  the  crystal  structure  and  the  average  simulation  structure  show  FMP 
strongly  interacting  with  polar  groups  of  the  active-site  cavity.  Most  important  are 
several  backbone  carbonyls  that  point  toward  the  center  of  positive  charge  density  of 
the  formycin  ring  and  make  hydrogen  bonds  to  its  protons.  In  addition,  both  structures 
show  N-6  donating  a  hydrogen  bond  to  the  carbonyls  of  Val  81  and  Gly  121,  and  a 
hydrogen  bond  between  N-7  and  the  carbonyl  of  Gly  121.  N-1  is  shown  to  accept  a 
hydrogen  bond  from  the  backbone  amide  group  of  Val  81. 

Tne  average  simulation  structure  is  in  further  agreement  with  the  crystal 
structure  showing  FMP  interacting  with  charged  residues  of  the  active  site,  princ'paliy 
Glu  177,  Arg  180,  an:.  Glu  208.  Residues  177  and  180  are  located  near  the  C- 
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terminus  of  an  a*helix  which  contains  a  distinct  bend  allowing  the  side  chains  of  both 


residues  to  reach  the  molecular  surface  of  the  active-site  cavity  From  the  simulation 
structure  the  two  residues  form  an  ion-pair  interaction  similar  to  the  interactions 
observed  in  the  native  enzyme  crystal  structure  [20,21,22]  and  the  full  MD  simulation 
structure  [23].  The  guanidinium  group  of  Arg  180  donates  a  strong  hydrogen  bond  to 
N-3  of  the  formycin  ring.  From  the  simulation  structure  the  carboxylate  of  Glu  177  and 
the  backbone  carbonyl  of  Glu  208  make  strong  electrostatic  interactions  with  the 
ribose  0-2'  and  0-3',  respectively.  From  the  X-ray  structure  the  ribose  0-2'  appears  to 

be  weakly  hydrogen  bonded  to  both  glutamates. 

TABLE  1 

COMPARISON  OF  THE  INTRAMOLECULAR  AND  INTERACTION  ENERGIES  OF 
FMP  AND  APG  SUBSTRATES  FOR  THE  X-RAY  STRUCTURE  AND 
SIMULATION  MODELS 


Energy 

FMP 

ApG 

(kcal/mol) 

X-ray* 

Simulation  t> 

X-ray* 

Simulation  b 

Intramolecular 

angle 

51.6 

47.6  (0.5) 

28.8 

47.7  (2.9) 

dihedral 

11.9 

26.0  (3.5) 

25.6 

29.3(1.9) 

bond 

44.8 

19.1  (5.7) 

5.6 

23.5  (2.9) 

X-term 

0.5 

2  6  (0.2) 

0.4 

0.6  (1.2) 

electrostatic 

-346.3 

-340.9  (2.0) 

-465.5 

-465.8(1.7) 

van  der  Waals 

7.3 

7.2  (2.7) 

14.2 

14.5(3.1) 

total 

-230.2 

-238.4  (6.2) 

-390.9 

-350.2  (11.9) 

Interaction 

electrostatic 

-68.1 

-173.9(1.8) 

-26.0 

-139.9  (8.5) 

van  der  Waals 

-13.8 

-15.3(1.9) 

-15.2 

-15.6(0.2) 

total 

-81.9 

-189  2(0.1) 

-41.2 

-155.5  (8.7) 

•  Values  for  the  X-ray  structure  calculated  using  partial  atomic  charges  of  the 

corresponding  simulation  model. 

b  Simulation  values  listed  were  averaged  over  200  ps  for  FMP  and  300  ps  for  ApG. 
the  values  in  parentheses  indicate  the  pertinent  standard  deviation. 
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The  simulation  structure  and  crystal  structure  show  contributions  to  the  bind.ng 
of  the  ligand  by  the  sioe  chain  of  lie  172  and  amide  groups  of  Asn  122,  Trp  211,  and 
Gly  212.  lie  172  packs  against  the  formycin  ring  and  makes  a  favorable  contribution 
through  van  der  Waals  interactions.  The  amide  backbone  of  residues  122,  211,  212 
exert  their  effect  on  the  binding  of  FMP  primarily  through  long-range  electrostatic 
interactions  with  N-7  for  residue  122,  and  both  the  ribose  0-3'  and  phosphate  group 
for  residues  21 1  and  212. 


Figure  1.  Stereoview  overlay  showing  a  superposition  of  the  average  simulation 
structure  and  X-ray  crystal  structure  for  the  binding  RTA-FMP  complex. 

Structural  differences  between  the  X-ray  crystal  structure  and  the  average 
simulation  structure  for  the  molecular  model  system  are  derived  primarily  from 

interactions  between  the  phosphate  group  and  protein  atoms.  From  tre  X-ra, 
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structure  the  phosphate  group  projects  out  cf  the  binding  pocket  into  the  solvent  and 
makes  a  long-range  (4.7  A)  electrostatic  interaction  with  the  positively  charged 
guanidinium  group  of  Arg  258.  The  average  simulation  structure  shows  the  phosphate 
group  at  a  lower  position  as  a  result  of  strong  ion-pair  interactions  with  Arg  258  and 
long-range  interactions  with  Arg  180.  Nonbonoed  interaction  energy  for  residue  258 
corresponding  to  the  simulation  structure  is  approximately  two-fold  greater  than  the 
crystal  structure. 

The  average  simulation  structure  shows  the  phosphate  group  strongly 
interacting  with  eight  water  molecules,  three  of  which  are  bridge  forming  waters  with 
side  chains  of  Tyr  80,  Gly  212,  and  Arg  258.  Because  of  strong  electrostatic 
interactions  of  the  phosphate  group  with  protein  atoms,  the  simulation  structure 
suggests  a  less  restrictive  solvent  accessible  surface  area  for  the  ribose  than  observed 
in  crystal  structure  showing  two  possible  hydrogen  bonds  with  water  molecules,  one 
bridging  Gly  212  and  the  phosphate  group.  The  simulation  also  clearly  shows  the 
desolvation  of  several  key  active-site  residues  upon  binding  cf  FMP.  A  comparison  of 
the  RTA-FMP  simulation  structure  with  the  average  MD  simulation  structure  of  the 
native  enzyme  shows  the  loss  of  approximately  five  water  molecules  bound  to 
residues  Glu  177,  Arg  180,  and  Glu  208.  These  waters  are  thought  to  be  critical 
nucleophiles  in  the  development  of  a  transition-state  like  binding  complex  of  the 
ribosome  [8,20,23]. 

Examination  of  the  structures  of  FMP  in  the  RTA  active  site  during  the  200-ps 
molecular-dynamics  simulation  shows  noticeable  thermal  fluctuations  in  the 
phosphate  backbone  and  formycin  ring  while  the  sugar  remained  relatively  rigid  along 
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the  nbose  0-4'  backbone  and  somewhat  flexible  in  the  hydroxyl  groups. 


TABLE  2 

COMPARISON  OF  X-RAY  AND  COMPUTED  TORSION 
ANGLES  FOR  FMP  IN  THE  RTA  ACTIVE  SITE 


Torsion 

angle* 

X-ray 

Simulation 

a 

166.2 

169.2  (6.3)b 

P 

-173.9 

88.8  (8.1) 

7 

-70.1 

-110.6  (5.4) 

5 

75.2 

86.1  (5.0) 

00 

44.1 

42.1  (3.7) 

0i 

-40.4 

-33.0  (4.8) 

02 

21.5 

12.5  (6.4) 

03 

6.5 

13.8(5.9) 

04 

-31.2 

-35.2  (4.4) 

X 

128.1 

113.2  (5.9) 

■  Employing  IUPAC  nomenclature  for  backbone  torsion 
angles,  ribose  endocyclic  torsions,  and  glycosidic  torsion. 
b  Values  in  parentheses  indicate  the  pertinent  standard 
deviation. 


Table  2  lists  a  comparison  of  the  torsion  angles  describing  the  conformation  of 
FMP  obtained  from  the  X-ray  structure  and  the  average  simulation  structure.  The 
calculated  backbone  torsion  angles  (3  and  y  show  deviations  from  the  crystal  structure 

of  about  94°  and  41°,  respectively,  and  can  be  attributed  to  the  significant  interactions 
between  protein  atoms  and  the  phosphate  group  of  the  simulation  model.  The 
simulation  structure  shows  exocyclic  and  endocyclic  torsions  corresponding  to  a  3'- 
endo  pucker  in  agreement  with  the  conformation  observed  in  the  crystal  structure. 
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Calculated  endocyclic  torsions  appear  to  be  driven  by  interactions  of  the  ribose  ring 
with  glutamates  1 77  and  208.  From  the  crystal  structure  the  formycin  ring  shows  a 

high  syn  configuration  with  a  torsion  %  of  about  128°  and  the  simulation  model  shows 
X  having  a  value  near  1 13°  ±  6°. 

Several  structural  motifs  constructed  from  the  original  FMP  structure  were 
investigated  using  free-energy  simulation  methods  for  the  dianion  and  monocnion 
models.  For  the  2-amino  formycin  5'-phosphate  analog  transformation,  the  free 

energies  calculated  for  the  enzyme-bound  ligand  (AAtind)  predict  the  analog  to  interact 
more  favorably  with  the  enzyme  in  both  dianion  and  monoanion  simulation  models. 
However,  the  relative  binding  free  energies  (AAA)  predict  the  2-amino  analog  to  bind 

less  strongly  to  RTA  than  does  FMP  itself.  Selectivity  reflects  solvation  effects,  as  it  is 
less  difficult  to  desolvate  FMP  than  its  2-amino  derivative.  The  lack  of  significant 
sensitivity  in  the  relative  free-energy  change  on  the  ionization  state  of  the  phosphate 

group  suggests  a  dominate  van  der  Waals  component  contribution  to  AAA. 

Of  the  free-energy  simulations  for  the  2-hydroxyl  formycin  5'-phosphate  analog,  the 
monoanion  model  shows  a  greater  relative  binding  affinity  (AAA  =-1.8  kcal/mol)  for  the 

active  site  of  RTA  than  the  original  FMP  substrate.  Analysis  of  the  average  end-point 
structure  shows  several  changes  in  the  structural  interactions  between  the  ligand  and 
residues  Tyr  80,  Glu  177,  and  Arg  180.  Most  important,  the  carboxylate  of  177 
interacts  favorably  with  the  hydroxyl  group.  The  formycin  ring  of  the  hydroxyl  analog 


16 


appears  to  be  a  high  syn  conformation  (x=  150°)  and  makes  several  strong 

interactions  found  in  the  FMP  simulation  structure,  most  notably  interactions  with  the 
polar  backbone  of  Val  81  and  Gly  121  as  well  as  packing  interactions  with  Tyr  80  and 
Tyr  123.  Finally,  decomposition  of  AAA  shows  a  large  electrostatic  free-energy 
interaction  component  dominating  the  improved  binding  of  the  ligand. 

As  for  the  2'-amino  formycin  5'-phosphate  analog,  neither  of  the  two  simulation 

models  predict  a  relative  increase  in  binding  affinity  for  the  enzyme.  Both  AAbina  and 

AAgoiy  show  a  decrease  in  favorable  interactions  with  the  enzyme  and  solvent, 
respectively. 

Molecular-dynamics  simulation  of  the  structural  interactions  governing  the 
binding  of  ApG  with  the  active  of  RTA  was  carried  out  for  300  ps.  Summarized  in 
Table  1  are  the  averaged  intramolecular  and  interaction  energies  for  ApG  in  the  ricin  A- 
chain  active  site  computed  from  the  simulation  and  the  corresponding  values 
determined  from  the  crystal  structure.  Energies  are  partitioned  into  covalent, 
electrostatic,  and  van  der  Waals  terms  in  accord  with  the  form  of  the  AMBER  potential 
function.  The  energy  contributions  that  show  significant  deviations  are  the  angle  and 
bond  strain  terms  as  well  as  the  electrostatic  interaction  binding  term.  Figure  2  shows 
a  superposition  of  the  average  simulation  structure  and  X*ray  crystal  structure.  The 
root-mean-square  (RMS)  deviation  for  the  enzyme-substrate  complex  is  roughly  0.66 
A  and  a  fit  of  the  substrate  structures  yields  1.21  A.  From  the  crystal  structure,  the 
adenine  ring  is  conformationally  stacked  between  the  rings  of  tyrosines  80  and  123 
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Both  residues  provide  efficient  packing  interactions  stabilizing  adenine  binding.  In 
addition,  both  tyrosines  stabilize  the  binding  of  the  guanosine  ribose  by  making 
several  long-range  electrostatic  interactions  and  van  der  Waals  contacts.  Comparable 
hydrophobic  contacts  are  found  in  the  average  simulation  structure,  however,  the 
substrate  conformation  is  shifted  by  approximately  1  to  3  A  with  the  largest  Cartesian 
coordinate  differences  occurring  in  the  adenine  and  guanine  rings.  This  shift  is  driven 
in  part  by  a  strong  hydrogen  bond  between  the  amide  backbone  of  Tyr  123  and  N-7  of 
adenine.  The  crystal  structure  shows  no  strong  enzyme  ion-pair  interactions  with  the 
imidazole  ring;  the  amide  backbone  of  Tyr  123  is  at  a  distance  of  3.62  A  and  functions 
to  secure  both  hydroxyl  groups  on  the  guanosine  ribose.  The  width  of  the  adenine 
binding  «*lot,  as  measured  from  the  distance  between  the  centroids  of  the  tyrosine 
rings,  is  about  9.93  A  in  the  simulation  structure  versus  8.86  A  in  the  crystal  structure, 
and  9.15  A  in  the  substrate-free  A-chain  crystal  structure  [20].  Centroids  measured 
from  the  pyrimidine  ring  to  the  tyrosine  rings  in  the  simulation  structure  are  5.49  A  and 
5.25  A  for  tyrosines  80  and  123,  respectively,  and  3.77  A  and  6.15  A,  respectively,  in 
the  crystal  structure.  The  packing  arrangement  of  the  tyrosines  with  the  adenine  ring 
in  the  simulation  model  exhibits  greater  edge-to-face  interactions  characteristic  of 
aromatic-aromatic  interactions  in  proteins  described  by  Burley  and  Petsko  [24],  A 
comparison  of  rotameric  states  for  residues  80  and  123  in  the  average  simulation 
structure  and  crystal  structure  finds  noticeable  variances.  For  Tyr  80,  Xi  =  (-145°, 

simulation;  -153°,  X-ray)  and  X2  =  (-132°,  simulation;  82°,  X-ray),  and  for  Tyr  123,  =  (- 

155J,  simulation;  -81°,  X-ray)  and  X2=  (-64\  simulation;  59°,  X-ray).  It  is  important  to 
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note  that  tne  side-chain  of  Tyr  80  undergoes  several  structural  transitions,  as 
discussed  further  below.  In  spite  of  differences  in  orientation,  the  net  contribution  of 
the  two  tyrosines  to  the  total  van  der  Waals  interaction  energies  binding  the  adenine 
ring  is  calculated  to  be  about  -4  kcal/mol  in  both  the  simulation  structure  and  the 
crystal  structure. 


Figure  2.Stereoview  showing  a  superposition  of  the  average  simulation  structure  and 
X-ray  crystal  structure  for  the  binding  RTA-ApG  complex. 

In  addition  to  the  hydrogen  bond  between  N-7  and  Tyr  123,  the  average 
simulation  structure  shows  the  adenine  ring  partially  “solvated"  by  :on-pair  interactions 
between  the  protons  at  N-6  and  the  backbone  carbonyls  of  Val  81  and  G!y  121.  The 
crystal  structure  shows  similar  interactions  for  N-6  plus  a  strong  hydrogen  bond 
between  N-1  and  the  amide  group  of  Val  81  (2.49  A  versus  4.64  A  in  the  simulation). 
The  significance  of  these  substrate-protein  interactions  with  the  polar  backbone  of 
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residues  81  and  121.  is  that  they  conceivably  provide  the  hydrogen-bonding  pattern 
that  specifically  recognizes  bases  in  the  order  ApG  versus  GpA,  which  has  been 
shown  net  bind  to  the  A-chain  [8].  Examination  of  the  dihedral  angles  d  and  \jr  for  Val 

81  shows  the  simulation  structure  located  in  the  same  region  of  conformational  space 
as  the  crystal  structure  placing  the  residue  at  (-45°, 76°)  compared  with  (-56°, 88°).  In 
contrast,  Gly  121  in  the  X-ray  crystal  structure  exhibits  dihedral  angles  of  (-62°, -97°), 
while  the  simulation  samples  the  conformation  (60°,-119°).  From  the  free  enzyme 
crystal  structure  [20],  Gly  121  adopts  the  conformation  (-82°,-6°)  and  contains  a 
hydrogen  bond  between  the  carbonyl  and  the  hydroxyl  of  Tyr  80.  This  hydrogen  bond 
is  replaced  by  substrate  interactions  with  the  net  effect  of  Val  81  and  Gly  121 
stabilizing  binding  by  -3  keal/mol  (-2  keal/mol  electrostatic)  in  the  simulation  structure 
and  -4  keal/mol  (-2  keal/mol  electrostatic)  in  the  crystal  structure. 

Substrate  binding  in  the  simulation  structure  and  the  crystal  structure  is 
significantly  stabilized  by  the  guanidinium  group  of  Arg  180  which  forms  several  ion- 
pair  interactions,  one  of  which  is  a  key  hydrogen  bond  with  N-3  cf  the  adenine  ring. 
The  simulation  structure  predicts  a  distance  of  about  2.85  A  between  the  guanidinium 
group  and  N-3  compared  with  2.55  A  observed  in  the  crystal  structure.  Mutagenic 
studies  have  implicated  Arg  180  as  the  principal  catalytic  residue.  The  conversion  of 
Arg  180  to  a  lysine  showed  minor  effect  on  enzyme  activity,  while  the  conversion  to  a 
histidine  resulted  in  a  1000-fold  reduction  in  activity  [25].  From  a  mechanistic 
viewpoint,  Arg  180  together  with  a  water  molecule  is  thought  to  be  critical  for 
protonating  adenine  at  N-3  in  facilitating  bond  breakage  during  deadenylation  of  the 
rRNA  substrate  [8.20].  Both  the  average  simulation  structure  and  the  crystal  structure 
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show  tha  substrate-bound  active  site  lacking  a  water  positioned  near  Arg  180.  *nough 
the  refined  crystal  structures  of  RTA  and  the  recombinant  A-chain  exhibits  c  'rapped 
water  bound  to  the  guanidinium  group  and  Glu  177  [20.22].  The  crystallographic 
thermal  factor  for  this  water  is  roughly  8  A2[20]  and  reflects  its  thermodynamic  stability 
in  the  structure  of  the  active  site.  The  loss  in  solvent  contribution  to  stah-lizing 
substrate  interactions  with  Arg  180  is  compensated  by  the  formation  of  ion-pair 
interactions  between  the  positively  charged  guanidinium  group  and  the  ?  denosine 
ribose  0-2'  together  with  favorable  long-range  electrostatic  interactions  with  the 

guanosine  and  phosphate  group.  There  appears  to  be  little  difference  in  the 
orientation  of  the  Arg  side  chain  in  the  free  and  substrate-bound  structures  (simulation 
and  X-ray),  both  showing  a  (— .t)  rotameric  state.  Likewise,  there  is  little  difference  in 
solvent  accessible  surface  area.  The  simulation  predicts  for  Arg  180  a  total  interaction 
energy  with  ApG  of  roughly  -38  kcal/mol  and  the  crystal  structure  yields  -31  kcal/mol. 
Energetic  differences  are  derived  primarily  from  the  phosphoryl  and  guanosine  ribose 
interactions.  Decomposition  of  the  Arg  180  structural  interaction  energies  in  terms  of 
substrate  functional  groups  for  the  simulation  structure  finds  -1 1  kcal/mol  interaction 
with  the  adenine  ring,  -16  kcal/mol  with  the  adenosine  ribose.  -9  kcal/mol  with  the 
phosphoryt  group,  and  *2  kcal/mol  for  the  guanosine  ribose  structure. 

The  highly  conserved  residue  Glu  177  has  been  proposed  to  stabilize  the 
transition-state  oxycarbonium  ion  on  adenosine  in  the  ribosome  substrate  [8). 
Analysis  of  the  conformation  of  ApG  from  both  the  crystal  structure  and  the  average 
s:"i'i!a"on  structe'e  shews  the  adepine  nng  in  a  high  syn  cenrieu-at  cn 
further  below)  lacking  the  required  orientation  to  mimic  cation  stabilization  on  the 
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nbose  by  Glu  177.  From  the  crystal  structure  the  carooxylate  shows  a  hydrogen  bora 
with  the  adenosine  0-2'.  while  the  simulation  structure  shows  no  direct  ion-pa 

interactions  with  the  substrate,  but  contains  a  strong  hydrogen  bond  with  Arg  180  as 
was  found  in  the  free  enzyme  (20).  Although  Glu  177  is  an  important  residue  in  the 
enzyme  activity  of  the  A-chain  [26].  an  examination  of  the  interaction  energies  of  both 
structures  suggests  a  net  destabilizing  contribution  to  the  binding  of  ApG  by  the 
carboxytate  primarily  through  poor  interactions  with  adenosine  ribose  in  the  simulation 
structure,  and  the  adenine  ring  and  guanosine  ribose  in  the  crystal  structure. 
Presumably  the  carboxylate  plays  a  key  role  in  anchoring  Arg  180  for  substrate 
binding  and  catalysis  as  well  as  maintaining  the  structural  integrity  of  the  active  site. 
Recent  molecular-dynamics  simulations  of  point  mutations  Glu  177  to  Gin  and  Ala  m 
the  native  ricin  A-chain  indicates  significant  destabilizing  tree-energy  electrostatic 
interactions  arising  from  the  guanidinium  group  of  Arg  180  in  the  mutant  folded 
proteins  [27]. 

Other  key  residues  making  favorable  substrate  binding  interactions  in  the 
average  simulation  structure  include  Glu  127,  Arg  134.  and  Arg  213.  Of  these 
residues,  only  Arg  134  contributes  lo  substrate  binding  in  the  c-ystal  structure.  The 
carboxylate  of  Glu  127  lies  near  the  C  terminus  of  an  a-helix  and  appears  to  stabilize 

the  binding  of  the  guanosine  ribose  through  long-range  electrostatic  interactions  with 
the  hydroxyl  groups  (calculated  to  be  about  -5  kcal/mol).  The  side  cham  of  the  highly 
conserved  residue  Arg  134  is  located  within  a  large  solvent  exposed  cod  running 

along  the  active-site  deft  ana  makes  several  lavoraDie  eectrcstat.c  interactions  w:  • 
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the  substrate;  the  adenine  ring  (-8  kcal/mol.  simulation-  -2  kcal/mol.  X-ray),  both  ribose 
structures  (-20  kcal/mol.  simulation;  -12  kcal/mol.  X-ray),  and  the  phosphate  group  (-12 
kcaLmci  in  both  simulation  and  X-ray).  Similar  strong  electrostatic  interactions  are 
observed  tor  Arg  213  in  the  simulation  structure  w.th  the  phosphate  backbone  (-11 
kcal/mol)  and  the  guanosme  bose  (-12  kcal/mol).  both  of  which  are  lacking  in  the 
crystal  structure.  The  guamdinn-.r  group  of  Arg  213  protects  out  from  an  a-helix  with  a 

predicted  (-.♦)  rotamer  conformation,  while  the  crystal  structure  shows  a  (-.t) 
conformation  with  the  residue  exposed  more  towards  the  solvent  but  less  than  in  the 
free  A-cham  crystal  structure  (20)  Collectively  these  predicted  protein-substrate 
interactions  *or  arginines  134  and  213  in  the  simulation  model  of  the  RTA-ApG 
complex  contribute  -63  kcal/mol  toward  stabilizing  substrate  binding  and  may 
conceivably  mimic  certain  individual  recognition  points  for  the  RNA  stem  and  loop 
residues  flanking  the  adenine  ring,  although  Arg  213  is  poorly  conserved  in  the  N- 
glycosidase  toxin  family.  Monzmgo  and  Robertus  [8]  have  constructed  several 
theoretical  models  binding  the  hexanucleotide  sequence  CGAGAG  with  'he  active  site 
of  ncm  A-cham  employing  the  NMR  defined  •  •  CGAAA  ••  hairpin  sequence  determined 
by  Heus  and  Pardi  (28).  From  their  models,  both  arginines  appear  to  display  ion-pair 
interactions  which  are  charactenstic  of  the  type  of  structural  elements  found  in  the  ApG- 
bound  RTA  simulation  structure.  In  one  of  their  models.  Arg  134  ion-pairs  with  the 
pnosphate  backbones  of  G2  and  A3  (the  target  of  deadenyiation).  and  Arg  213  pairs 
with  the  phosphate  of  Cl.  Further  stabilizing  interactions  are  found  for  Arg  258  with 

The  simuiat.on  structure  «s  w  agreement  with  the  crystal  structure  showing  two 
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hydrogen  bonds  between  the  guanosine  moiety  and  residues  Asp  124  and  Asn  209. 
The  amide  backbone  of  Asp  124  makes  an  ion-pair  interaction  with  the  ribose  0-2' 

(2.03  A,  simulation;  2.02  A,  X-ray)  and  the  side-chain  amide  nitrogen  of  conserved  Asn 
209  donates  a  hydrogen  bond  to  a  phosphoryl  oxygen  (1.81  A,  simulation;  2.61  A,  X- 
ray).  The  simulation  structure  predicts,  however,  a  net  12  kcal/mol  greater  binding 
stabilization  for  both  residues  than  observed  in  the  crystal  structure.  Additional  protein 
interactions  with  the  substrate  are  formed  by  the  conserved  residue  He  1 72  and  the 
invariant  residue  Trp  211.  The  side  chains  of  172  and  211  make  favorable 
hydrophobic  interactions  with  the  adenine  and  adenosine  ribose,  respectively.  The 
indole  ring  lies  parallel  and  provides  addition  anchoring  support  for  the  side  chain  of 
Arg  180.  Chemical  modifications  of  Asn  209  to  a  Ser  and  Trp  21 1  to  a  Phe  have  been 
shown  to  produce  little  effect  on  the  depurination  of  the  ribosome  [26]  suggesting 
involvement  in  base  recognition  or  related  functions. 

Binding  differences  between  the  average  simulation  structure  and  the  crystal 
structure  can  be  further  characterized  from  analyzing  the  substrate  solvent  accessible 
surface  area  and  molecular  volume.  The  crystal  structure  shows  the  substrate  more 
exposed  to  solvent  with  a  25  A2  greater  surface  area  (employing  a  1.4  A  radius 
spherical  probe)  and  a  300  A3  greater  molecular  volume.  Taken  together,  the 
contraction  of  the  average  simulation  substrate  structure  at  protein  surface  is 
consistent  with  stronger  binding  interactions.  The  simulation  model  shows  the 
substrate  interacting  with  five  water  molecules.  The  diffusion  coefficient  of  the  bound 
waters  calculated  using  the  Einstein  relation  is  roughly  0.14  As  ps-1,  compared  with 
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0.70  As  ps-'  for  all  the  waters  in  the  simulation  system.  Two  of  the  bound  waters  3re 
bridge-forming  waters  with  guanosine  and  protein  side  chains;  one  bridging  0-5'  and 
Arg  213,  and  the  second  bridging  H03'  with  residues  Asp  124  and  Arg  134.  One  of 
the  waters  appears  to  be  bound  by  the  substrate  in  a  four-sided  negatively  charged 
density  region  composed  of  N-3,  0-4',  0-5'  of  guanosine  and  0-4'  of  adenosine.  The 

remaining  waters  are  hydrogen  bonded  to  protons  at  N-1  and  N-2  of  the  guanosine. 

As  shown  in  Table  I,  the  X-ray  crystal  structure  and  the  average  simulation 
structure  are  comparable  in  total  intramolecular  electrostatic  energies.  The  crystal 
structure  shows  a  strong  intramolecular  hydrogen  bond  between  H0-3'  of  the 
guanosine  ribose  and  a  phosphoryl  oxygen.  Weaker  hydrogen  bonds  may  be  formed 
between  HO-3'  and  either  0-4'  of  adenosine  or  0-5'  of  guanosine.  The  simulation 

structure  contains  two  hydrogen  bonds;  a  strong  hydrogen  bond  between  the  hydroxyl 
group  of  the  adenosine  ribose  and  a  phosphoryl  oxygen  plus  a  second  weaker  bond 
between  0-5'  of  adenosine  and  a  proton  at  N-2  of  the  guanine  ring.  The  simulation 

showed  little  variation  in  the  bond  common  to  the  hydroxyl  group  (1.7  ±  0.1  A), 
whereas  the  adenosine-guanine  hydrogen  bond  appears  to  be  iess  stiff  (2.5  ±  0.3  A). 
For  comparison  purposes  the  intramolecular  distances  between  C-5'  and  N-9  for  both 
nucleotides  as  a  function  simulation  time  have  been  computed.  The  simulation 
oredicts  averaged  distances  of  3.9  ±  0.4  A  and  4  8  ±  0.1  A  for  adenosine  and 
guanosine,  respectively,  versus  3.9  A  and  4.3  A,  respectively,  for  the  crystal  structure. 
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The  simulation  shows  a  transition  of  about  t  A  in  the  adenosine  C*5'-N-9  distance 

which  may  be  facilitated  in  part  by  a  short-lived  conformational  transition  in  the  nbose 
(discussed  further  below).  The  average  RMS  fluctuation  for  ApQ  calculated  for  the 
heavy-atoms  is  about  0.48  A  compared  with  1.33  A  estimated  from  the  crystallographic 
thermal  factors  of  the  refined  structure  [8].  The  guanosine  base  and  ribose  showed  the 
largest  thermal  fluctuations  consistent  with  lower  bindu  j  interaction  energies. 

A  comparison  of  the  torsion  and  pseudorotation  angles  describing  the  substrate 
conformation  obtained  from  the  X-ray  crystal  structure  and  the  average  simulation 
structure  is  listed  in  Table  3.  The  backbone  torsional  angles  a  and  (3  predicted  by  the 

simulation  show  deviations  from  the  crystal  structure  of  about  30°  and  94°, 
respectively,  and  can  be  attributed  to  differences  in  electrostatic  interactions  between 
the  phosphate  group  and  arginines  134  and  213  as  well  as  intramolecular  hydrogen 
bonds.  The  calculated  dihedral  angles  yand  5  for  adenosine  are  in  good  accord  with 

the  crystal  structure,  whereas  the  simulation  predicts  a  different  conformation  for 
guanosine.  Similarly,  the  predicted  endocyclic  torsions  of  the  adenosine  ribose 
compares  well  against  the  crystal  structure,  unlike  most  guanosine  values.  Both  the 
average  simulation  structure  and  the  crystal  structure  exhibit  syn  base  configurations 
for  the  adenine  and  guanine  rings,  however  the  conformations  differ  by  roughly  57°. 

Analysis  of  the  crystal  structure  shews  C3 '-exo  and  C3 '-endo  ribose  pucker 
configurations  for  adenosine  and  guanosine,  respectively.  The  simulation  model 
indicates  that  the  adenosine  ribose  is  in  a  dynamic  equiiibr.-.n  ccr:  sting  of  C3'-e*o 
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and  C2 '-endo  configurations.  Short-lived  dihedral  transitions  to  a  Cl'-exo 


configuration  occurred  during  the  simulation.  The  greatest  fluctuations  were  found  in 
dihedrals  0i.  02.  and  0a  for  adenosine.  The  adenine  ring  exhibited  a  glycosidic  torsion 
fluctuation  of  about  50°  at  the  latter  transition,  but  remained  in  syn  configuration.  No 
transitions  occurred  for  guanosine,  which  showed  a  C4'-exo  pucker  configuration 
throughout  the  simulation. 


TABLE  3 

COMPARISON  OF  X-RAY  AND  COMPUTED  TORSION  ANGLES  FOR  ApG  IN  THE 
RTA  ACTIVE  SITE 


Torsion 

angle 

X-ray 

Simulation 

A 

G 

A 

G 

a 

-158.0 

-187.7  (8.7) 

P 

-81.1 

-175.0(13.3) 

7 

•168.4 

-47.6 

-175.3(13.3) 

41.8(10.2) 

8 

144.2 

81.4 

144.1  (8.8) 

65.4  (5.8) 

00 

-33.1 

43.3 

-32.6  (4.1) 

23.1  (5.4) 

0i 

33.5 

-42.5 

30.0(19.9) 

-40.1  (4.6) 

33.5 

24.5 

-22.5(10.6) 

42.0(5.1) 

03 

-1.4 

4.5 

3.9(10.1) 

-27.2  (6.4) 

@4 

22.5 

-31.5 

20.5  (6.1) 

1.2  (6.5) 

X 

94.2 

15.5 

150.8  (10.3) 

73.4  (9.3) 

p 

196.6 

13.1 

198.2 

57.2 

•  Values  in  parentheses  indicate  t..e  pertinent  standard  deviation. 


To  further  characterize  the  conformational  dynamics  of  ApG  in  the  active  site  of 
the  A-chain,  side-chain  and  backbone  dihedral  angles  of  several  key  residues  were 
examined  as  a  function  of  simulation  time.  The  side  chains  of  Tyr  80  and  Arg  ISO 
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snowed  significant  dihedral  angle  changes  which  were  greater  than  the  fluctuations  of 
the  angle.  Of  these  residues,  Tyr  80  clearly  showed  angle  changes  which  directly 
influence  and/or  accommodate  the  conformational  changes  in  the  bound  substrate. 
Structurally,  Tyr  80  is  the  initial  residue  in  an  antiparailel  8-sheet  following  a  reverse 
turn,  and  as  shown  above,  provides  a  thermodynamically  stable,  hydrophobic 
environment  for  the  substrate.  Relative  to  the  X-ray  crystal  structure  of  the  free  enzyme 
[20],  Tyr  80  is  the  only  residue  which  shows  any  appreciable  conformational  change 
upon  binding  ApG.  The  crystallographic  dihedral  angles  (X1.X2)  are  displaced  from 

(-151°.-86°)  to  (-153°, 80°).  Initially  from  the  trajectories,  xi  and  *2  are  located  at 

(-86°, -130°).  At  25  ps  xi  is  displaced  and  oscillates  around  a  slightly  different  value 

from  the  initial  equilibrated  conformation,  yet  the  side  chain  remains  within  the  same 
rotameric  state.  This  subtle  change  appears  to  be  in  response  to  the  60°  or  greater 
increase  in  61  of  the  adenosine  ribose.  The  trajectories  showed  Tyr  80  undement 

large  angle  excursions  at  roughly  50  ps  through  100  ps,  at  165  ps.  210  ps,  and  275  ps. 
Fluctuations  starting  at  90  ps  results  in  a  transient  packing  defect  leading  to  the 
transition  at  100  ps  with  the  potential  minima  located  at  (-175°,-116°).  The  transition  at 
165  ps  correlates  with  the  short-lived  Cl'-exo  repuckering  in  the  adenosine  ribose 

coupled  with  the  fluctuations  in  the  adenine  glycosidic  torsion.  There  are  two 
additional  transitions  at  about  210  ps  and  275  ps.  The  former  transition  is  due  to  the 
replacement  of  a  hydrogen  bond  between  the  backbone  carbonyl  of  Tyr  80  and  the 

backbone  amide  of  Asp  75  with  a  bond  betweer,  the  carbonyl  and  the  am:ne  group  of 
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the  adenine  ring.  The  average  dynamic  structure  shows  a  CO-HN  (adentne)  distance 
of  4.62  A,  as  compared  with  4.35  A  from  the  X-ray  crystal  structure.  The  latter 
transition  in  corresponds  to  fluctuations  in  the  hydrogen-bonding  pattern  of  the  side- 
chain  atoms  for  Asp  75  and  Asn  78. 

The  simulation  predicts  several  structural  transitions  in  the  backbone  angles  of 
Tyr  80  which  are  correlated  with  the  calculated  fluctuations  in  adenosine  outlined 
above.  From  0  to  50  ps.  (<>,v)  has  oscillations  near  (-52°,  1 14°),  which  represents  the 
observed  crystallographic  conformation.  This  state  is  followed  by  large  fluctuations  up 
to  100  ps  leading  to  the  second  state  where  (p.vy)  has  oscillations  at  (19°, 81°),  and  is 

interrupted  by  the  short-lived  transition  state  at  165  ps  with  oscillations  near 
(-128°,  130°).  The  fourth  state  has  oscillations  at  (53°, -43°),  remaining  throughout  the 
simulation.  No  similar  transitions  were  found  in  the  simulation  for  Tyr  123  or  any  other 
ac*ive-site  residue  nor  were  there  any  comparable  transitions  found  in  the  MD 
simulation  of  the  free  enzyme  [23].  Thus,  suggesting  Tyr  80  may  function  as  the  key 
determinant  of  substrate  modulation  and  orientation  in  the  active  site  once  recognition 
has  occurred,  principally  by  arginines  134,  180,  and  213.  Relative  to  the  free  enzyme 
simulations  [23],  ’he  magnitudes  of  the  side-chain  and  main-chain  positional  RMS 
fluctuations  for  Tyr  80  in  the  bound  simulation  are  significantly  altered.  There  also 
appears  to  be  substrate  effects  on  the  dynamics  of  Asp  75  plus  reverse-turn  residues 
76-79.  The  main-chain  and  side-chain  fluctuations  of  these  residues  change  by  as 
much  as  0.2-0.4  A  between  the  free  and  bound  simulations.  The  free  enzyme 
structure  compr.ses  cf  a  reverse-turn  hydrogen  bond  between  HN  of  Asn  76  ana  ire 
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carbonyl  of  Asp  75  together  with  a  (3-sheet  bond  between  the  carbonyl  of  Tyr  80  and 
HN  of  Asp  75.  As  shown  above,  this  latter  ion-pair  interaction  switches  during  the 
bound  simulation  to  pair  with  the  adenine  ring.  The  average  simulation  structure 
exhibits  several  ion-pair  interactions  between  Asn  78  and  various  residues  within  the 
reverse-turn  and  (3-sheets;  the  backbone  carbonyl  of  Asn  78  pairs  with  the  amide 
group  of  Tyr  80;  the  carboxylate  of  Asp  75  pairs  with  the  nitrogen  side-chain  of  Asn  78 
plus  the  amide  group  of  Thr  77;  and  the  amide  backbone  groups  of  Asn  78  and  Ala  79 
pair  with  the  carbonyl  of  Asp  75.  The  X-ray  crystal  structure  cf  the  bound-state  shows 
similar  hydrogen  bonds  with  the  notable  exception  of  forming  an  alternative  ion-pair 
interaction  between  the  amide  group  of  Tyr  80  and  the  side-chain  oxygen  of  Asn  78. 
These  additional  reverse-turn  interactions  observed  in  the  bound-state  structures 
appear  to  be  stabilized  by  the  binding  of  ApQ.  It  is  reasonable  that  this  structural 
element  corresponding  to  residues  75-80  could  play  an  important  entropic  role  in 
binding  substrates. 

The  dihedrals  of  Arg  180  showed  large  fluctuations  near  135  ps,  but  remained 
in  a  (-,t)  rotameric  state.  Fluctuations  were  of  a  jump-and-return  nature  and  appear  to 
be  due  primarily  to  shifts  in  hydrophobic  interactions  with  Glu  177  and  Trp  211.  There 
were  no  corresponding  fluctuations  in  the  backbone  angles  for  Arg  180.  For  arginines 
134  and  213,  the  trajectories  exhibited  no  fluctuations  in  134  and  only  minor  side- 
chain  displacements  in  213,  which  were  correlated  with  the  adenosine  phosphate 
backbone  y  displacements.  Comparisons  with  the  free  enzyme  simulation,  arginines 

134,  180,  and  213  showed  decreased  positional  fluctuations  in  the  bound-state 
simulation.  The  loss  of  flexibility  in  these  key  residues  in  the  presence  of  ApG  is 
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readily  understood  in  terms  of  specific  enzyme-substrate  interactions  outlined  above. 

It  is  of  interest  to  compare  the  binding  of  ApG  with  that  of  FMP.  As  mentioned  in  the 
Introduction,  analysis  by  X-ray  crystallography  suggests  that  ricin  A-chain  has  a 
binding  preference  for  FMP  [8].  An  examination  of  the  binding  motifs  of  the  adenine 
and  pyroazolopyrimidine  rings  in  the  two  substrates  shows  similarity  in  packing  and 
ion-pair  formation,  however  there  are  several  positional  and  conformational  shifts. 
These  adjustments  correspond  primarily  to  the  pyroazolopyrimidine  ring  able  to  form 
an  extra  hydrogen  bond  to  the  polar  backbone  of  the  active  site  coupled  with  fewer 
electrostatic  and  orientation  constraints  of  the  phosphate  backbone  in  FMP.  A 
comparison  of  the  total  binding  interaction  energies  calculated  from  the  MD  simulation 
of  the  RTA-FMP  complex  with  the  simulation  of  the  RTA-ApG  complex  indicates  a 
difference  of  about  40  kcal/moi,  favoring  FMP  binding.  Both  simulations  predict 
significantly  greater  binding  interactions  than  observed  in  the  crystal  structures. 
Differences  in  structural  interactions  between  the  two  substrates  are  primarily 
electrostatic  in  origin  as  hydrophobic  contacts  are  roughly  equal.  The  average 
simulation  structure  of  the  ApG  complex  contains  less  efficient  electrostatic  binding 
interactions,  principally  those  of  the  dinucleotide  phosphoryl  group  and  adenosine. 
When  combined  with  the  guanosine  ribose,  the  two  ribose  structures  of  ApG 
approximate  the  FMP  ribose-protein  interactions.  Key  residues  which  make  significant 
contributions  to  the  binding  of  FMP  but  not  ApG  include  Asp  124.  Glu  208,  Gly  212. 
and  Arg  258.  The  first  two  residues  make  favorable  electrostatic  interactions  with  the 
ribose  group;  the  backbone  carbonyl  of  208  forms  a  strong  hydrogen  bond.  The 
amide  backbone  of  Gly  212  and  Trp  211  together  with  the  guamdmium  group  of  Arg 
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253  (which  substitutes  for  Arg  213  in  the  ApG  binding  mode)  provides  partial 
neutralization  of  the  negative  charge  of  the  phosphoryl  group.  As  in  the  binding  of 
ApG,  Arg  180  makes  significant  stabilizing  interactions  securing  FMP  in  the  active  site 
and  forms  a  strong  ion-pair  interaction  with  N-3.  The  carboxylate  of  Glu  177  makes 
favorable  binding  electrostatic  interactions  with  the  ribose  (unlike  that  of  ApG).  but 
poor  interactions  with  the  pyroazolopyrimidine  ring. 

Building  on  the  average  simulation  structures  foi  substrates  FMP  and  ApG 
bound  in  the  RTA  active  site,  design  studies  are  being  carried  out  aimed  at 
determining  new  ligands  with  improved  binding.  Current  efforts  are  directed  at 
appending  new  substituents  onto  the  existing  FMP  and  ApG  structures.  Shown  below 
(structure  III)  and  in  Figure  3  is  an  example  of  a  proposed  ligand  structure  constructed 
from  molecular-dynamics  simulations  of  the  adenine-like  ring  in  FMP  bound  to  the 
active  site  of  RTA. 


Ill 
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CONCLUSIONS 


Ricin  A -chain  is  a  cytotoxic  protein  that  attacks  ribosomes  by  hydrolyzing  the  adenine 
ring  from  a  specific  adenosine,  thereby  inhibiting  protein  synthesis.  The  substrate 
binding  complexes  of  formycin  5 '-monophosphate  and  adenyl-3',5'-guanosine  in  the 

ricin  A-chain  active  site  were  investigated  by  molecular-dynamics  simulations 
immersed  in  solvent  water.  The  studies  undertaken  showed  the  average  simulation 
structures  of  the  substrate-bound  enzyme  to  be  in  good  accord  with  the  observed  X- 
ray  crystal  structures  in  reproducing  overall  binding  modes.  However,  there  are 
significant  differences  in  the  location  and  binding  of  various  substrate  functional 
groups.  For  FMP,  the  simulation  predicts  significant  differences  in  the  location  and 
binding  of  the  phosphate  group.  For  ApG,  the  calculated  average  structure  of  the 
bound  substrate  is  shifted  1  to  3  A  from  the  X-ray  structure,  while  retaining  a  similar 
binding  motif  for  adenosine.  The  simulation  predicts  greater  protein-ApG  binding 
interactions,  particularly  among  the  guanosine  ribose  and  phosphoryl  function  groups. 
Free-energy  simulation  methods  have  been  employed  to  explore  several  structural 
motifs  of  FMP  which  would  have  a  greater  binding  affinity  for  the  active  site.  It  is  shown 
that  ricin  A-chain  has  a  preference  for  FMP  over  analogs  2-  and  2 '-amino  formycin  5  - 
phosphate  and  2-hydroxyl  formycin  5 ‘-phosphate.  Using  the  binding  motif  of  the 
adenine  ring  from  the  average  simulation  structures,  several  substituents  have  been 
appended  to  the  base  with  the  removal  of  the  ribose  and  phosphate  group  leading  to 
the  design  of  new  ligands  for  ricin.  These  potential  antidotes  are  being  further 
evaluated  by  molecular-dynamics  simulations  to  determine  the  relative  binding 
affinities. 
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